Excitonic effects and optical properties of passivated CdSe clusters 
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We calculate the optical properties of a series of passivated non-stoichiometric CdSe clusters using 
two first-principles approaches: time-dependent density functional theory within the local density 
approximation, and by solving the Bethe-Salpeter equation for optical excitation with the GW 
approximation for the self-energy. We analyze the character of optical excitations leading to the 
first low-energy peak in the absorption cross-section of these clusters. Within time-dependent density 
functional theory, we find that the lowest-energy excitation is mostly a single-level to single-level 
transition. In contrast, many-body methods predict a strong mixture of several different transitions, 
which is a signature of exciton effects. We also find that the majority of the clusters have a series 
of dark transitions before the first bright transition. This may explain the long radiative lifetimes 
observed experimentally. 



Experimental advances in the synthesis of semiconduc- 
tor clusters have stimulated considerable theoretical ef- 
fort to understand the optical and electronic properties 
of these systems. Semiconductor clusters often exhibit 
strong size-dependent effects, which are not yet fully 
understood. As an intermediate system between sin- 
gle atoms and bulk materials, semiconductor clusters are 
also of intrinsic theoretical interest. Clusters of II- VI el- 
ements, such as CdSe, have attracted considerable atten- 
tion in recent years owing to their potential technological 
applications in various devices such as solar cells, lasers 
and biological imaging tools, among others 0. Great 
effort has gone into fabricating and characterizing size- 
controlled samples. This is particularly challenging be- 
cause clusters of different sizes have similar stoichiometry 
and are synthetized by similar reactions in which the tem- 
perature, the solvent and the ratio and concentration of 
precursors have to be carefully controlled |2| . Theoretical 
calculations are difficult because of the inherent complex- 
ity of accurate theories, as we can infer from the limited 
number of theoretical articles in the literature |3|, Li 13, |f| . 

The use of pseudopotentials and density functional 
theory (DFT) have been very successful in determining 
the ground-state properties of both bare 3] and passi- 
vated 5] CdSe clusters. However, DFT is a ground- 
state theory, and it has serious shortcommings in provid- 
ing a quantitative description of optical and electronic 
excitations Understanding spectroscopic experi- 

ments requires the computation of excited state prop- 
erties, which present a greater challenge than ground- 
state calculations. The problem is addressed by both 
time-dependent density functional theory (TDDFT) and 
GW/Bethe-Salpeter (GW/BSE) methods. TDDFT is 
simpler to implement, but a good general approxima- 
tion for the exchange-correlation functional is still lack- 
ing . The local-density approximation within TDDFT 
(TDLDA) has been found to give accurate results for 



some finite systems such as sodium clusters j| . In oth- 
ers, it gives at best a qualitative picture 0,00- On the 
other hand, GW/BSE has been shown to be very accu- 
rate in bulk materials , albeit more computationally 
demanding. Until recently, first-principles GW/BSE cal- 
culations have been done only for very small clusters, 
containing no more than 35 atoms [y, [ll| . 

In this letter we present a detailed comparison of 
TDLDA and GW/BSE calculations of the optical prop- 
erties of passivated CdSe clusters. By examining the 
character of the transition leading to the first peak in 
the absorption spectra, we investigate the importance 
of many-body effects, fully accounted for in GW/BSE 
but absent in TDLDA. Calculations of energy band gaps 
and absorption spectra of bare CdSe clusters have been 
done using TDDFT within the local-density approxima- 
tion (TDLDA) 3] . Energy band gaps of passivated CdSe 
clusters have also been calculated within the TDDFT 
framework 5]. GW/BSE has been used for calculations 
of optical properties of molecular systems, silicon clus- 
ters, and in crystals 00 El E3> but not ni " V or n ~ 
VI clusters yet. Experimental data on CdSe clusters is 
readily available |2| which makes these systems ideal for 
a comparative analysis. 

We studied a series of five clusters: Cd4Se6, CdgSe^, 
CdioSei6, Cdi7Se28 and Cd.^Se so. While larger CdSe 
clusters are found to be spherical , these smaller clus- 
ters are of pyramidal shape (see Fig. [U). Two of the 
clusters studied have zincblende structures (Cd4Se6 and 
CdioSei6), while the rest are of the wurtzite type. All 
clusters were passivated by fictitious hydrogen atoms of 
charge 1.5e (attached to surface Cd atoms) and 0.5e (at- 
tached to surface Se atoms) [l4j , in order to simulate the 
effect of surfactants on the surface of the real clusters . 

The initial geometry of the clusters was constructed 
based on the X-ray data 2] and then relaxed [01 EH- 
Tetrahedral symmetry was conserved during the relax- 
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FIG. 1: Geometry of the clusters studied: Cd4Se6, CdsSei3, 
CdioSei6, Cdi7Se28 and Cd32Se5o- Cd atoms are dark yellow, 
Se atoms are light yellow, and fictitious H atoms are small 
gray. 



The matrix elements for Q are given by 



(2) 

where Xuw = nw — ™kr are the difference between the oc- 
cupation numbers, and Ttuoikr = £kr — is the difference 
between the eigenvalues of the single-particle states. K is 
the coupling matrix which describes the linear response 
of the system |2lj . 

Electron-phonon coupling and temperature depen- 
dence effects are included a posteriori by broadening the 
absorption spectra with a normalized Gaussian function 
with fixed dispersion of 0.1 eV. 

In the GW/BSE method, the many-body expression 
for the polarizability II is related to the electron-hole 
correlation function L by: 



ation. The Cd-Se bonds in the final structures are com- 
parable to the experimental measurements 2] , but about 
5% shorter than the Cd-Se bond in bulk cadmium se- 
lenide. Calculations were done using norm-conserving 
pseudopotentials [Tjj constructed within the local density 
approximation of density functional theory. Pseudopo- 
tentials for cadmium and selenium have scalar relativis- 
ts effects included and the interaction between the Cd 
Ad orbital and valence orbitals is accounted for by a non- 
linear core correction in the pseudopotential Spin- 
orbit and semicore effects beyond the non-linear core cor- 
rection are ignored. In order to estimate how important 
these effects are, we computed the band gap of bulk CdSe 
in the wurtzite structure. The gap obtained within the 
GW approximation is 1.8 eV, which compares well with 
the spin-orbit averaged gap obtained in experiment: 1.97 
eV |l9j. The measured spin-orbit splitting, A, at the 
valence band maximum at the T point is A = 0A3eV 
. Including spin-orbit splitting would change the cal- 
culated energy gap by E gap = y ~ 0.1 5eV. The Kohn- 
Sham equations are solved on a real s pac e grid using a 
higher-order finite difference method 20]. All calcula- 
tions were done within a spherical boundary of radius at 
least 6 a.u. from the outermost passivating atoms. A grid 
spacing of 0.3 a.u. was used for LDA calculations, while 
a grid spacing of 0.6 a.u. was used for the calculation of 
optical properties as described below. 

In TDLDA 0, [2^ , the optical response is evaluated 
as a first-order perturbation in the electron density due 
to an external potential. The excitation energies Q n are 
obtained from a solution of the eigenvalue equation 0, 
El: 



QF n = Q 2 n F n (1) 



n(l,2) = -*L(l,2;l+ 2+) (3) 
where L satisfies the Bethe-Salpeter equation 0,03]: 

L(1,2;3,4) = G(1,4)G(2,3) + 
J d(5678)G(l, 5)G(6, 3)K(5, 7; 6, 8)L(8, 2; 7, 4). (4) 

Solving the above equation, we obtain optical excita- 
tions of the electronic system. G is the electron Green's 
function, and the kernel operator K describes interac- 
tions between the excited electron and the hole left be- 
hind in the electron sea. The electron self-energy is calcu- 
lated within the GW approximation [2^. As in TDLDA, 
the absorption spectra is broadened by normalized Gaus- 
sian functions with fixed dispersion. In both TDLDA and 
GW/BSE, the optical gap is defined as the energy of the 
first transition with measurable oscillatory strength (the 
first allowed transition) [22^ . 

Figure shows optical gaps as a function of CdSe clus- 
ter size. For all but the smallest cluster, our calculations 
p show a trend very similar to that found in experi- 
ment. The dependence of the gap on cluster size is very 
strong. It is interesting to note that while TDLDA cal- 
culations underestimate the gap by ~0.5 eV to ~1.5 eV, 
GW/BSE overestimates the gap by less than 0.6 eV. The 
discrepancy between experimental data and TDLDA cal- 
culations increases as the size of the cluster increases, 
but the opposite is observed for GW/BSE calculations. 
Based on the analysis made for bulk CdSe, we expect 
that neglected spin-orbit and semicore effects at pseu- 
dopotential level are responsible for a residual discrep- 
ancy between theory and experiment of a few tenths of 
electron-volt. 

There is a large discrepancy between the theoretical 
calculations and the experimental measurements for the 
optical gap of the smallest cluster, Cd4Se6- For the larger 
clusters, all TDLDA results are below experiment, while 
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FIG. 2: Experimental and calculated optical gaps. The LDA 
gap is simply the difference between Kohn-Sham eigenvalues. 



BSE results are above. This is not the case for Cd4Se6, as 
both TDLDA and BSE predict larger optical gaps than 
measured experimentally. Given the reduced number of 
atoms in this cluster, it is not clear whether the organic 
ligands on the surface or some other mechanism might 
be responsible for this discrepancy. 

Calculated absorption cross-sections (normalized by 
the total number of cadmium and selenium atoms in 
each cluster) are shown in Fig. 03 Both TDLDA and 
BSE have a well-defined first peak at low energy, with 
a second peak separated from the first by ~1 eV. We 
have analyzed the character of the excitation leading to 
the first peak observed in the absorption cross-section. 
As shown in Tabled the first TDLDA excitation in all 
clusters is mostly a result of single-level to single-level 
transitions. For the three smallest clusters, the dominant 
transition is from the HOMO (highest occupied molecu- 
lar orbital, which is triple-degenerate for most clusters, 
without spin-orbit splitting) to the LUMO (lowest un- 
occupied molecular orbital, which is non-degenerate in 
all cases). This is not so for the two largest clusters as 
they have a series of dark transitions (transitions with 
negligible oscillatory strength) before the first bright (al- 
lowed) transition. The first optically allowed transition 
for Cdi7Se28 involves transitions from the third level be- 
low the HOMO to the LUMO. For Cd 3 2Se 5 o the first al- 
lowed transition is from the first level below the HOMO 
to the LUMO. 

Experimentally, bulk CdSe in the wurtzite and 
zincblende structures is a direct gap semiconductor [2^ . 
Our own LDA calculations agree with these experimen- 
tal results. Clusters, however, do not always behave 
in the same way. This can be attributed to confine- 
ment effects and the geometry of each cluster 0, Eij . 
Cd4Se6, CdgSeis and CdioSei6 have dipole allowed tran- 
sitions between the triple-degenerate HOMO and the 
non-degenerate LUMO. But the HOMO of Cdi 7 Se 28 is 
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FIG. 3: Absorption cross section calculated within TDLDA 
(dashed curve) and GW/BSE (solid curve). 



non-degenerate because of a change in the ordering of 
the energy levels (the HOMO-1 is triple-degenerate), 
and the HOMO-LUMO transition becomes dipole for- 
bidden as a result of selection rules. Cd32Se5o has a 
dipole allowed HOMO-LUMO transition, but the oscil- 
lator strength is small because of little overlap between 
HOMO and LUMO wavefunctions. 

The character of the GW/BSE excitation leading to 
the first peak in the absorption spectra is very different 
from the one predicted by TDLDA. For the smallest clus- 
ter the excitation is still dominated by a single-level to 
single-level (HOMO -> LUMO) transition. For the rest 
of the clusters, however, the excitation is the result of 
a strong mixture of different transitions. There are two 
sources of mixing in GW/BSE: at the GW level, mixing 
occurs because of fact that LDA wavefunctions are not 



TABLE I: Energy of the first allowed transition, E; presence 
of dark transitions (those with negligible oscillatory strength) 
before the first allowed transition; and percentage of the low- 
est energy peak in the absorption cross section that is due to 
the single- level to single- level transition indicated in the text. 
No entry in the percentage column indicates that the transi- 
tion is strongly mixed, i.e, the largest component contributes 
less than 25%. 
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4.77 


no 
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5.97 


yes 


91 


CdgSeis 
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no 
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no 
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3.21 
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96 


4.45 
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Cdi7Se28 


2.38 


yes 
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3.81 


yes 
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1.87 
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3.65 
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identical to quasiparticle wavefunctions |g, la l2i| 5 at the 
BSE level, the electron-hole kernel is stronger and more 
non-local than the TDLDA kernel. We find that the mix- 
ing in TDLDA is one order of magnitude smaller than 
the mixing in GW/BSE. This is a result of stronger cou- 
pling matrix elements at the GW/BSE level compared to 
TDLDA. We quantify the mixing as 



FG02-89ER45391 and DE-FG02-03ER15491. Calcula- 
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Institute, the Texas Advanced Computing Center and at 
the National Energy Research Scientific Computing Cen- 
ter (NERSC). 



(5) 



where H is the effective Hamiltonian matrix in either 
TDLDA or GW/BSE methods. M = corresponds to 
a situation where there is no mixing between different 
single electron transitions. We find that the mixing M 
within GW/BSE is of the order of 10 -2 while within 
TDLDA it is of the order of 10" 3 . In both cases the 
mixing is non-zero but still much smaller than one. As a 
result of the stronger mixing in GW/BSE, the excitation 
cannot be associated with a single electron-hole transi- 
tion. A similar behavior has been observed in small sili- 
con clusters and in bulk semiconductors 0, , and it is 
a signature of excitonic effects. 

The effects of mixing also explain the observed diver- 
gence in the experimental and TDLDA curves in Fig. 
Since the energy levels in the smaller clusters are more 
separated, the mixing effect is not as large and TDLDA 
calculations are more accurate than for larger clusters 
where the energy levels are closer together and more mix- 
ing can occur. 

In the GW/BSE calculations we find that there are 
dark transitions for all clusters but CdgSe^. Van Driel 
et al. 4] have recently shown that measured rates of 
emission are completely determined by radiative decay 
and that the occupation of dark excitonic states consid- 
erably attenuates spontaneous emission. The presence of 
dark transitions in our GW/BSE calculations then may 
explain in part the long radiative lifetimes (~l-10/is) ob- 
served experimentally. 

In conclusion, we have calculated the optical proper- 
ties of a series of small CdSe clusters using two different 
approaches: TDLDA and GW/BSE. We find that the 
two methods lead to a very different character for the 
lowest energy excitation. In TDLDA, the excitation is 
dominated by a single-level to single-level transition. In 
GW/BSE, however, the excitation is the result of strong 
mixing between different transitions. We interpret this 
as due to exciton effects. Our calculations also show 
that most clusters have a series of dark transitions before 
the first bright transition, which attenuates spontaneous 
emission and may explain the long radiative lifetimes of 
these clusters. 
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